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Abstract. - The interaction of a quantum system with a bath, usually referred to as dissi- 
pation, can be controlled if one can establish quantum interference between the system-bath 
interaction and a coupling of the system to an external control field. This is demonstrated for 
the example of the spin-boson model in the strong coupling limit for the system-bath inter- 
action. It is shown that driving and trapping of the spin system leads to an optimum control 
problem which is nonlinear in the external control field. Using an indirect optimization strategy 
introducing a Lagrange-type adjoint state, we show that the spin system can be trapped in 
otherwise unstable quantum states and that it can be driven from a given initial state to a 
specified target state with high fidelity. 



Introduction. - Over a last a few years, a number of interesting schemes have been pro- 
posed to eliminate the undesirable effects of decoherence in open quantum systems, including 
decoherence free subspaces [1,2], quantum error correction codes [3-5], quantum feedback 
[6] and mechanisms based on unitary "bang-bang" pulses and their generalization, quantum 
dynamical decoupling [7,8]. The key ingredient of dynamical decoupling is the continuous 
disturbance of the system, which suppresses the system-environment interaction. It has been 
shown that, in the bang-bang control schemes, the decoherence of the system is effectively 
suppressed if the pulse rate is much higher than the decoherence rate due to the system- 
environment interaction [7]. 

The staring point of the decoupling techniques is the observation that even though one 
does not have access to the large number of uncontrollable degrees of freedom of the envi- 
ronment, it is still possible to interfere with its dynamics by inducing interactions into the 
sub-system which drive it so fast that the environment cannot follow [8]. Alternatively, if 
one can establish a suitable coupling to the system by means of an external control, one can 
establish quantum interference with the system-bath. In a simple minded model for a dissi- 
pative quantum system, where the interference between the system-bath and system-control 
interaction is ignored or is irrelevant only limited control can be achieved [9]. The situation 
changes dramatically when interference between the system-environment and system-control 
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interaction can be used to control the effective system-environment coupling [7,8,10-19]. The 
degree and nature of quantum interference constructive or destructive can be controlled by 
adjustment of the control field, known as coherent control. 

In this paper we apply the concept of coherent control to steer a dissipative quantum 
system to the spin boson model, in which a quantum two-level system (qubit) is modelled by 
a spin, the environmental heat bath by quantum oscillators and the spin subjected to external 
control field is coupled to each bath oscillator independently [19 -21]. Achieving decoherence 
control for this model is formulated using optimal control which is mathematically a problem 
of functional optimizations under constraints in form of differential equations [9,22]. The 
two-level system coupled to a bath provides an adequate model of such diverse phenomena as 
electron transfer reaction [17], elect ron-phonon interaction in point defects [23] and quantum 
dots [24], interacting many-body systems [25], magnetic molecules [26] and bath assisted 
cooling of spins [12]. 

Bloch-Redfield formalism. - Consider a physical system S embedded in a dissipative 
environment B, also referred to as the heat bath, and interacting with a time-dependent 
classical external field i.e., the "control". The total Hamiltonian iJtot = Hs{t) + Hb + i?int 
is the sum of the system Hamiltonian Hs{t), the bath Hamiltonian Hb and their interaction 
Hint, which is responsible for decoherence. Note that the operator Hs{t) contains a time- 
dependent external field to control the quantum evolution of the system. We suppose that 
the system-environment interaction Hamiltonian is bilinear i?int = X]q ® B,^ where 
and Ba are Hermitian operators of the system and the environment, respectively. 

The basic assumptions underlying the derivation of the equation of motion for the reduced 
density matrix p{t) = Tyb {ptot(i)} , are that (i) the initial factorization ansatz; we assume 
that at time t = the bath B is in thermal equilibrium and uncorrelated with the system S 
(ptot(O) = p(0)(8)pB, Feynman- Vernon approximation), (ii) weak system-bath interaction limit 
in which the second-order perturbation theory is applicable (ptot(i) = p{t) ® PB + 0{Hnt), 
Born approximation) (iii) the relaxation time tb of the heat bath is much shorter than the 
time scale tr over which the state of the system varies appreciably {tb <S tr, justifying the 
Markov approximation). From the Liouville-von Neumann equation ihptot — [-^tot , Ptot] for 
the total density operator and after performing the above mentioned approximations, one 
obtains the master equation for the reduced density matrix in Bloch-Redfield form 



where the first term represents the unitary part of the dynamics generated by the system 
Hamiltonian Hs{t) and the second term accounts for dissipative effects of the coupling of the 
system to the environment. The Redfield relaxation tensor Rijki{t) is given by [27] 




(1) 
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where the time-dependent rates F^.j.^ {t) are evaluated through the following expressions: 
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Fig. 1 - Population transfer from the ground state ^(0) = 1 to the maximally mixed state z{tf) = 0. 
(a) shows the control field vs. time, (b) and (c), respectively, show F^^ and Fo vs. time, (d) shows 
the relative population z — (az) vs. time. Wf — 1, Wr — 0, a = 0.2, eo — A, iVc = 20A, /3 = 1/hA 
and tf = 100/A. The guess parameters chosen are {A, il, 0) = (A, A/10, 7r/3) and those computed 
read {Aopt, fiopt, <l>opt) = (5.838 A, 0.211 A, 28.173 rd). A is an arbitrary unit of frequency. 



Here C/s(i, t') = T |cxp —j^J^, (It Hs{t) | is the propagator of the coherent system dynamics 

while {Ba{t)B/}{0)) = Tyb {B{t)B{0)pB} is the environment correlation function with pB = 
exp{—/3HB)/Z, the canonical ensemble of the bath at the inverse temperature (3 = I/UbT. 

Because of the applied control field, the transition rates defined by Wji{t) = ^tj jii^) + 
^ij ji{t) , in the secular approximation which we suppose also valid in the driven case, become 
time dependent [15,16]. The field influence on both the unitary and dissipative contributions to 
the time evolution of the physical system makes possible an external control of dissipation. In 
particular, the correlation between the control field and the dissipation leads to the destruction 
of the detailed balance limt^oo Wij{t)/Wji{t) ^ exp{—f3Ei)/ex.p{—l3Ej) where Ei are the 
energy eigenvalues of the undriven physical system. So, the steady state can be far from 
equilibrium in the driven case. The influence of the control field on the relaxation tensor via 
Us{t,t') is a direct consequence of quantum interference between the system-bath interaction 
and the coupling of the system to the external field. 

Driven spin boson model. ~ The Hamiltonian of the driven spin boson model where the 
two-level system is bilinearly coupled to an ensemble of harmonic oscillators is given by [19,21] 

E = -^A<T^ - ^ (^0 + ^(*)) + mujfxf] + (a, qo/2) ^ . (5) 

where aa with a = x,y,z are Pauli spin matrices; HA is the tunneling splitting, Hsq is an 
energy bias and fi€{t) is its modulation by a time-dependent external control field; and the 
heat bath is represented by a set of harmonic oscillators of mass mi, angular frequency Wj, 
momentum pi and position coordinate Xj. The oscillators are coupled independently to the 
spin coordinate with strength measured by the set {ci} while go measures the distance between 
the left and right potential wells. The coupling constants enter in the spectral density function 
of the environment defined by, J{ui) = ^ J2i T^TuF ^{^ — ^i) ■ 

In order to compute the Redfield tensor, it is necessary to determine the propagator of 
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the coherent system dynamics Us{t,T). An analytical expression for Us{t,T) is not trivial 
because the Hamiltonian of the physical system Hs{t) = —^Aa^ — f (^o + ctz is time- 
dependent and not diagonal. To get round this difficulty, we transform the Hamiltonian H by 
the unitary operator, i.e., polaron transformation [/ = with O = —^a^^ and fl = J2i 
where fij = [qoCi / hmioSf) pi [21]. The transformed Hamiltonian H' = UHU~^ takes the 
following form 

H' = - J {so + e{t)) a, + 1 V + m^oj^xA - ^hA (a+e*" + <7_e-*") , (6) 



where a± = (a^ iiay)/2. After invoking the polaron transformation the coherent propagator 
is trivial and reads Us[t, t') = cos [eo{t - t') /2 + f(t, f) /2] I + i sin [ea(t - f) /2 + /(t, t') /2] tj^ 
where the function J{t,t') = J^, dreir) captures the effects of the external control field. 

Rate master equation. - Let us consider the polaron transformed spin-boson Hamiltonian 
and derive the explicit form for the corresponding master equation for small A. Here, the 
system and the environment operators are Si = hAa+/2, S2 = S\ = fiA(T_/2 and Bi = e~*^, 
B2 = B\ = e*^, respectively. The Bloch-Redfield formalism leads to the following master 
equation for the diagonal elements of the the reduced density matrix (the populations) 

Pllit) = P22{t)Wi2it) - pii{t)W2l{t) , (7) 
P22{t) = Pll{t)W2l{t) - p22{t)Wi2it) , (8) 

with time dependent transition rates Wi2{t) = F^-^ -^glO + ^21 12 (^) ^^"^ W^2i(i) = 21 (^) + 
r^2 2i(*)- Up to the second order in A, the populations decouple from the non-diagonal terms. 
The master equation for the population difference z{t) = {(Jz)t = Pii(i) — P22{t) reads 

z{t) = -^zz{t)z{i) + To{t), (9) 

where T^^{t) = Wi2{t) - W2i{t) and Toit) = Wi2{t) + W2i{t) defined as 

T,,{t) = A^ f dre-«^(")cos[/(i,t-r)]cos[Qi(T)], (10) 

ro(t) = A2 / dre-«=(^)sin[/(i,t-T)]sin[Qi(T)]. (11) 
Jo 

The quantities Qi{t) and Q2{t) are the imaginary and the real part, respectively, of the func- 
tion $(r) = ll dio^ [(1 - coswr) coth(fex;/?/2) i sinwr)] where e-*(^) = (e"^(^)e-*"(°)) 
is the environment correlation function [28]. In the present work, the spectral density of the 
environment J{lo) = rjuje~'^/'^'' is assumed to be Ohmic with exponential cutoff. Here ry is 
a phenomenological friction coefficient. For an Ohmic bath and at low temperatures regime 
{hu}c(3 » 1) one can use the approximation [21] 

Qi{t) = 2q; arctan(a;cT), (12) 



Q2{t) = aln(l + w^r^)+2aln 



sinh(7rT/ ph) 
{nT/l3h) 



(13) 



where the dimensionless dissipation constant a = qo^/27r/i have been introduced. 
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Fig. 2 - Population transfer from the ground state z{0) = 1 to the target excited z{tf) = —1. (a) 
shows the control field vs. time, (b) and (c), respectively, show T^z and Fq vs. time, (d) shows the 
relative population z — (az) vs. time. Wf — 1, Wr — 0, a = 0.2, eo = A, Wc = 20A, f3 = 1/hA and 
tf = 100/A. The guess parameters chosen are {A,Q,(j)) = (A, A/10, 7r/3) and those computed read 
(v4opt,nopt,<Aopt) = (9.34A, -1.107 X 10~^A, 10.73 rd). A is an arbitrary unit of frequency. 



Quantum optimal control problem. - Let time t be in the interval [0,tf] for fixed tf. The 
evolution of the state variable z{t) governed by the the master equation depends not only on 
the initial state z{0) = Zi but also on the time-dependent control variable e{t). The task is 
now to find a control field that will steer the system from its initial state to a desired final 
state at specified time tf. Typically, it is possible to define a cost functional incorporating 
the objective. The goal of optimal control algorithms is to calculate a control field which can 
induce a specified system dynamics by minimizing this cost functional. Consider then the 
following quantum optimal control problem of minimizing the cost functional J(p) subject to 
the dynamical constraint, i.e, the master equation for z{t) = {(Tz)t 

min J(p) = ^{z{tf) - z,f + dt{z{t) - Zr{t)f , 

m = -TUt)z + ro{t), t€ [0, tf], (14) 
^(0) = Zi . 

The first term in the cost functional J(p) measures the deviation of the final state z{tf) from 
the desired state z^ . During the time interval [0,t/], there may also exist a desired state 
trajectory Zj.{t). This objective is incorporated in the second term. Here, it is assumed that 
< Wf < 1 and < Wr < 1 with Wf + Wr = 1. J(p) is the sum of the so-called final time 
cost functional and running cost functional. The vector p G is a set of parameters on 
which the control field depends, i.e., £{t) = e{t, p). Here, an optimal solution of this problem 
is characterised by first order optimality conditions in the form of the Pontryagin's Minimum 
Principle [22] . These conditions are formulated with the the help of a Hamilton function that 
has the following form in our problem; 

n{z. A, p) = ^{z{t) - zr{t)f + X{t) {-T,,{t)z + To{t)} . (15) 
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Fig. 3 - Trapping of the spin in unstable excited state, ^(0) = 1. (a) and (b), respectively, show 

Tzz and Fq vs. time, (c) shows the relative population z = (a^) vs. time, (d) shows the cost 
functional vs. the number of iteration. Wf = 0, Wr — 1, a — 0.2, eo = —A, u>c = 20A, /3 — 1/hA 
and tf = 50/A. The guess parameters chosen arc (A,Q,,(j)) — (12A, 5A,0) and those computed read 
{Aopt, fiopt, (t>opt) = (218.196A, 89.995A, -1.117 rd). A is an arbitrary unit of frequency. Note that 
Oopt is much higher than Wc- 



Pontryagin's minimum principle states that a necessary condition for {z,p) to be a solution 
of the above optimal control problem is the existence of an adjoint state A such that 



.Xt)z + To{t), 



z{Q) = Z„ X{tf) = Wr {z{tf) - Zd) , 

0=1^, 



(16) 



dpi ' 



l...Np and t€[0,tf] 



The minimum principle requires the solution of a set of complicated nonlinear algebraic 
equations, namely the optimality conditions dH/dpi = 0, i = l...Np which can only be 

solved in an iterative manner. The present optimal control problem is not singular because 

det (^ Qp.^. ^ 7^ 0, i,j ~ 1. . ■ Np, since Tzz{t) and Tait) depend nonlincarly on e(t,p). The 

optimality conditions dH/dpi = 0, i = 1 . . .N, gives the gradient for the cost functional J(p) 
with respect to control parameters pi 



dpi Jo 



dp., 



i dpi 



dt, 



l...N„ 



(17) 



In summary, the computation of the gradient towards an optimal solution requires the com- 
putation of z{t) forward in time and the adjoint state X{t) backward in time. 

Num,erical Results. - In the present work, we write the control field as monochromatic 
plane wave e(i,p) = e{t, A,^,(j)) = Acos{p.t + </>). In this case, an analytic expression can 
be found for the function /(t, t') = f{t, t', A, fl, (p) = J*, dTe{T, A, fl, </>) which greatly reduces 
the numerical effort for finding a solutions to this optimization problem. The task is then to 
find a set of three parameters namely the amplitude pi = A, the frequency p2 = ^ and the 
phase P3 = (f) such that the cost functional J(p) defined above is minimal. The optimization 
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is performed using the gradient method. More precisely, we used the subroutine DMNG of 
port library [29] implementing the quasi-Newton method as variant of gradient algorithms. 

Heating. - As a first example, we consider driving the system from a pure state z(0) = 1, 
corresponding to the temperature-zero ground state of the two-level system, into the mixed 
target state z{tf) = 0. Fig. 1(d) shows that in the absence of the control, the system inevitably 
relaxes to the thermal equilibrium state 

ro(eo) Wi2{(ii) - \V2i{co) . . . , , 

z^t = ^^-r^ = . X , „/ / X = tanh(fi/3eo/2) . (18) 
r^^(eo) W^i2(eo) + W2i(eo) 

Fig. 1 shows also guess and optimal harmonic control, obtained via the gradient method, the 
rate Tzz{t) and the inhomogeneous term To{t), as well as the time evolution of the relative 
population z{t). Owing to the control dependence of Tzz{t) and To{t) which reflects the 
periodic nature of the control field, Fig. 1(b) and 1(c), the objective posed can be achieved 
perfectly. 

Cooling. - As a second example, we consider a spin-flip in which (az) = 1 at t=0 
is changed to {(Tz) = —1 at target time tf with an optimal solution displayed in Fig. 2a. 
Here, the simplest optimum harmonic solution selected is a low-frequency field which simply 
changes the level sequence in the two-level system, at the same time, renormalising Tzz and 
Pp. Approaching target time tf, Tzz{t) approaches — ro(t) to provide z{tf) = —1, as shown 
in Figs. 2(b) and 2(c). In Fig. 2(d) we also show the instantaneous equilibrium value Zst(i) = 
tanh [hp (eo + Eopt (t)) /2] which corresponds to the ratio {W12 (t) - W21 (t)) / {W12 (t) + W21 (i)) 
= Toit)/^zz{t) relevant in the adiabatic limit, valid when \de{t)/dt\ <§; as discussed in the 
literature [19,26]. In the long-ime limit, the prediction from the adiabatic approximation is 
fulfilled by the numerically selected optimal electric field £opt{t) in Fig. 2(d). In the context 
of bath-assisted cooling of spins, our finding has certain analogies with the results of Ref. [12] 
where a version of the spin-boson model under influence of externally controlled pulses is stud- 
ied. Similar to our situation, an interference between the external fields and the spin-bath 
interaction create a mechanism that cool the spins much below the bath temperature. Starting 
from the initially infinite-temperature spin the authors cool it down to very low temperatures 
by increasing the spin's polarisation. 

Trapping . - As a third example, we consider trapping of the spin system for eq = —A 
in state z{t) = 1 for all time t G [0, i/] which now corresponds to the (unstable) excited state 
of the system. In isolated case, trapping is made possible by dynamic localization using a 
periodic control field to renormalize the coupling between the two levels from A to IS.Jq{A/VI) 
for a large frequency of the field; here Jo{x) is the zero order Bessel function [19,31]. Therefore, 
when A/Vlis such that Jo(j4/i7) — 0, the two-levels are decoupled and the tunneling between 
them is suppressed. In the presence of dissipation the trapping is destroyed, the undriven 
system relaxes to the equilibrium state, and the guess harmonic field with high frequency is 
not able to trap the system as shown in Fig. 3(c). For an optimum harmonic control field, 
we find that trapping becomes truly eS'ective when the control frequency exceeds the cut-off 
frequency which was set as uic = 20A. The optimum solution for the control parameters 
is (^opi.f^opt,</>opt) = (218.196 A, 89.996 A, -1.119 rad). As shown in Figs. 3(a) and 3(b), 
this high-frequency control renormalizes Tzz{t) and TQ(t) to values of the order 10~^A such 
that z{t) w for all time t & [0,t/]. This finding is in qualitative agreement with the bang- 
bang method used to study the independent spin-boson model [7]. As we mentioned in the 
introduction, this study showed that the control has to be switched at a rate greater than 
l/Wc in order to provide effective decoupling of the physical system from the environment. 
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Fig. 3(d) shows the cost functional versus the number of iteration. According to this figure, 
the cost functional is monotonically decreasing and its convergence to zero is reached after 15 
iterations. 

Conclusion. - In summary, we have formulated and solved an optimal control problem 
for the spin-Boson model to demonstrate feasibility of controlling the effective system-bath 
interaction by an external control field. An external control not only renormaliscs the spin 
Hamiltonian but also the effective coupling between system and environment. We have demon- 
strated that the effective system-bath interaction can be increased or decreased in controlled 
fashion. This mechanism can be used to optimize dynamic driving and trapping of the spin 
system. Results were presented for control of the relative population z of the spin system and 
a harmonic control field. This work has been extended to a control of all components of the 
Bloch vector simultaneously, as well as to general control field shapes. The results will be 
presented elsewhere. 

* * * 
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